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Calibration data of a wind tunnel sting balance was processed using a can- 
didate math model search algorithm that recommends an optimized regression 
model for the data analysis. During the calibration the normal force and the 
moment at the balance moment center were selected as independent calibration 
variables. The sting balance itself had two moment gages. Therefore, after ana- 
lyzing the connection between calibration loads and gage outputs, it was decided 
to choose the difference and the sum of the gage outputs as the two responses 
that best describe the behavior of the balance. The math model search algo- 
rithm was applied to these two responses. An optimized regression model was 
obtained for each response. Classical strain-gage balance load transformations 
and the equations of the deflection of a cantilever beam under load are used to 
show that the search algorithm’s two optimized regression models are supported 
by a theoretical analysis of the relationship between the applied calibration loads 
and the measured gage outputs. The analysis of the sting balance calibration 
data set is a rare example of a situation when terms of a regression model of 
a balance can directly be derived from first principles of physics. In addition, 
it is interesting to note that the search algorithm recommended the correct re- 
gression model term combinations using only a set of statistical quality metrics 
that were applied to the experimental data during the algorithm’s term selection 
process. 
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= lower integration interval bound 
= upper integration interval bound 
= coordinate of the center of the forward balance gage 
= coordinate of the center of the aft balance gage 

= distance between the center of the forward and the center of the aft balance gage 
= first auxiliary coordinate difference 
= second auxiliary coordinate difference 

= distance between model moment center and balance moment center 
= modulus of elasticity 

= force at the model moment center and at the balance moment center 
= force at the center of the forward balance gage 
= force at the center of the aft balance gage 

= distance between balance calibration fixture and balance reference axis 
= area moment of inertia 

= first regression coefficient of deflection angle equation at free end of sting 
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= second regression coefficient of deflection angle equation at free end of sting 
= corrected nominal distance between calibration force and balance moment center 
= nominal distance between the calibration force and the balance moment center 
= distance between calibration force and balance moment center after leveling 
= moment arm correction due to the force induced deflection of a beam 
= moment arm correction due to the moment induced deflection of a beam 
= moment at balance moment center 
= moment at model moment center 
= moment at the center of the forward balance gage 
= moment at the center of the aft balance gage 

= nominal moment computed with moment arm of unloaded sting balance 
= total number of terms of a regression model (includes intercept) 

= electrical output at the forward moment gage 
= electrical output at the aft moment gage 
= force location offset due to the calibration fixture design 
= coordinate of the sting balance reference axis 
= x-coordinate of a cantilever beam (zero at fixed end of beam) 

= regression coefficients 
= regression coefficients 

= constant describing relationship between electrical output and moment at balance gage 

= constant describing relationship between electrical output and moment at forward gage 

= constant describing relationship between electrical output and moment at aft gage 

= beam deflection due to a concentrated force 
= beam deflection due to a pure moment 

= total moment arm correction due to cantilever beam deflection 
= total moment arm correction due to elastic change of cantilever beam length 
= regression coefficients 

= theoretical estimate of regression coefficient e\ 

= regression coefficients 

= theoretical estimate of regression coefficient rji 
= theoretical estimate of regression coefficient 772 
= angular deflection 
= regression coefficients 
= coordinate of the balance moment center 

= number between the lower and upper bound of an integration interval 
= regression coefficients 
= standard deviation 

= coordinate of the balance calibration force 
= angle describing the location of the balance calibration force 

I. Introduction 


During the past 4 years a software package called BALFIT was developed for the Wind Tunnel Division 
at NASA Ames Research Center that is used for the regression analysis of wind tunnel strain-gage balance 
calibration data. The software uses an innovative candidate math model search algorithm in order to find 
an optimized regression model for the analysis of data. The algorithm is not limited to wind tunnel balance 
calibration data applications. It can also be applied to more general multivariate experimental data sets. 

In principle, the search algorithm tries to identify an optimized regression model, i.e., the so-called 
recommended math model, that has superior predictive capabilities. It is obtained after minimizing a 
search metric that is a measure of the regression model’s expected predictive capability. In addition, two 
constraints are applied that make it possible to test only those regression models during the search that 
meet strict statistical quality requirements. A detailed description of the search algorithm’s approach can 


2 

American Institute of Aeronautics and Astronautics 



be found in Refs. [1], [2], and [3]. 

So far, the search algorithm has been applied to a wide variety of wind tunnel balance calibration 
data sets. Traditional hand load data sets of multi-piece, single-piece, and floor balances were processed 
successfully (Refs. [4], [5]). In addition, large data sets from calibration machines were also analyzed (Refs. [6], 
[7]). A calibration data set of a wind tunnel sting balance was analyzed for the first time in 2008. This 
data set was courteously supplied to the authors by The Gulfstream Aerospace Corporation. During the 
analysis of the sting balance data set an interesting connection between the optimized regression model of 
the data and a first principles analysis of the relationship between the calibration loads and the gage outputs 
was observed. This connection is discussed in great detail in the present paper. At first, however, basic 
assumptions and elements of the regression model optimization algorithm are reviewed in order to provide 
a better understanding the regression model optimization process. 

II. Regression Model Optimization Algorithm 

Figure 1 shows basic elements of the most recent version of the candidate math model search algorithm. 
The goal of the algorithm is to find an optimized regression model, i.e., the so-called recommended math 
model, that best represents the given experimental data set. The regression model search (optimization) is 
performed in several steps. Initially, the user of the algorithm has to select a function class combination 
that appears to be most suitable for the global regression analysis of the experimental data. In the next 
step, the upper bound for the regression model search has to be identified. This upper bound is called the 
“permitted math model.” It is the largest possible math model for the given function class combination and 
data set that will lead to a non-singular solution of the global regression problem. This upper bound makes 
it possible to perform an automated regression model search as it is the largest regression model in the given 
context that wifi not lead to a software “crash” during the regression model search. 

How can the upper bound of the math model search space be identified? This question is answered 
in some detail in Ref. [3]. The basic approach of the identification process can be summarized in a few 
sentences. First, each math term of the regression model of the responses is considered to be a column 
vector that has a number of rows equal to the number of data points of the data set. A non-singular solution 
of the global regression problem can only exist if the column vectors representing the math terms are linearly 
independent. They have to be basis vectors of a vector space. 

At this point another question emerges. How can the linear independence of a set of vectors be tested? 
A numerical technique called Singular Value Decomposition (SVD) has to be applied “iteratively” for this 
purpose (see also the discussion in Ref. [3], p.5). What does “iteratively” mean in this context? Let 
us assume, for example, that a regression model of the responses consist of a total number of “n” terms 
(including the intercept term). Then, a total number of “n” column vectors can be defined that represent 
all possible math terms of a regression model of the data set for the given function class combination. Now, 
the first vector, e.g., the vector of ones representing the intercept term, is considered to be the initial basis 
vector of a vector space. In the next step, SVD is used to test if the second vector is linearly independent of 
the first vector. Only two possibilities exist. Possibility 1: The second vector is linearly independent of the 
first vector ==> a new vector space is defined using the first and second vector ==> the new vector space will 
be used to test the third vector. Possibility 2: The second vector depends linearly on the first vector the 
second vector is no longer considered to be a new basis vector of a vector space ==> the corresponding math 
term is removed from the regression model =>• the original vector space remains unchanged and will be used 
to test the third vector. The process is repeated for each vector. The vector space is increased, of course, by 
using the tested vector whenever this vector is linearly independent of the vector space that is defined by the 
previously found linearly independent set of vectors. Consequently, SVD has to be applied a total number 
of “n — 1” times in order to identify and remove all potentially linearly dependent vectors (math terms) that 
may be contained in the largest theoretically possible regression model. The remaining set of vectors (math 
terms) defines the final vector space for the given data set and function class combination that will always 
lead to a non-singular solution of the global regression problem. 

Theoretically, no order of the column vectors (math terms) has to be specified in order to test column 
vectors for linear dependency using the suggested iterative application of SVD. Lower order terms of the 
regression model may be tested before -or- after higher order terms. In a practical implementation, however, 
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testing lower order terms before higher order terms may provide a better understanding of the reason why 
a certain column vector (math term) was rejected during the iterative application of SVD. 

The search for the optimized regression model can begin after the upper bound of the math model 
search space has been identified by applying SVD to the selected function classes and data set. The search 
starts at the lower bound of the search space which, for wind tunnel strain-gage balance calibration data 
applications, is defined by the intercept and the principle linear load term that is connected with a specific 
balance gage output. For a sting balance, for example, the principle linear load term for the difference of 
the gage outputs is the normal force at the balance moment center (BMC). A metric, i.e., the standard 
deviation of the PRESS residuals of the responses is minimized during the search. This metric is used to 
compare the predictive capability of different regression models that are tested during the search (see also 
the flowchart in Fig. 1). A primary and secondary constraint are enforced during the search. The primary 
search constraint makes sure that a tested regression model only has statistically significant terms. The 
secondary search constraint helps avoid math models that have unwanted near-linear dependencies. In the 
end, only those regression models are considered during the search that fulfill both constraints (see Ref. [1] 
for more details about the search metric and the application of the two search constraints). 

The user has the option to apply the “hierarchy rule” to the optimized regression model after the com- 
pletion of the search. The “hierarchy rule” adds missing lower order terms to a regression model. The 
enforcement of the “hierarchy rule” is “optional” as its application will not necessarily lead to a regression 
model with better predictive capabilities (a detailed discussion of this viewpoint is given in Ref. [1]). Previ- 
ously, the search algorithm also allowed the user to enforce the “hierarchy rule” during the candidate math 
model search. This option is no longer supported as it appears to lead to suboptimal search results for 
certain types of data sets. 

In the next part of the paper the sting balance’s design and its calibration are discussed in more detail. 
Then, results of the regression model optimization of the calibration data will be presented. 

III. Wind Tunnel Sting Balance Design and Calibration 

The selected sting balance has two moment gages that are mounted forward and aft of the BMC. The 
drawing in Fig. 2 shows the location of the two moment gages relative to the BMC. The balance has an 
approximate total length of 28 [in] and a 30° bend angle at the midpoint. Therefore, the centerline of the 
test article is about 5 [in] below the centerline of the downstream portion of the balance. 

Figure 3a shows a sideview of the sting balance as it appeared during the calibration. The balance 
gages are covered using gage fairings. In addition, an Angle Measurement System (AMS) unit is attached 
to the test article mounting plate. The AMS unit was used to level the test article mounting plate during 
calibration. Figure 3b shows a more detailed view of the forward and aft gage after the removal of the gage 
fairings. 

The calibration of the sting balance was performed using weights that were applied in both normal and 
upside-down orientation of the balance. Figure 3c shows the balance, AMS unit, the calibration fixtures, 
and weights in upside-down orientation. This orientation was used to apply a positive normal force during 
the calibration. Four horizontal load points were used on the calibration equipment in order to vary the 
distance between the normal force and the BMC. Figure 4 shows all force and moment combinations that 
were applied during the calibration of the balance. Four distinct lines can be seen in Fig. 4. They are related 
to the four moment arms that were chosen for the calibration (16.52 [in], 18.52 [in], 20.52 [in], 25.92 [in]). 

The forward balance gage measures a raw electrical output called R\. The aft gage measures a raw 
electical output called R 2 . Both gages are moment gages. Consequently, the gage output R\ is approximately 
proportional to the moment Mi at the forward gage and the gage output R 2 is approximately proportional 
to the moment M 2 at the aft gage. Two principle combinations of regressors and responses are possible in 
order to perform a regression analysis of the calibration data. The two options depend on the balance load 
format that is selected for the regression analysis. The moments Mi and M 2 at the forward and aft gage 
may be used as independent variables for the regression analysis. In that case, the electrical outputs R\ and 
R 2 are good choices for the responses. However, it is also possible to select the force F and moment M at 
the BMC as independent variables. Then, the linear combinations R 2 — Ri and Ri + R 2 should be selected 
for the analysis as R 2 — Ri is nearly proportional to F and R\ + R 2 is nearly proportional to M. The two 
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options are summarized in Table 1 below. 


Table 1: Regressor and response combination options for sting balance calibration analysis. 



REGRESSORS 

RESPONSES 

OPTION 1 

Mi, M 2 

Ri, R 2 

OPTION 2 

F, M 

R 2 — Ri , Ri + R 2 


Why is the difference R 2 — RI nearly proportional to F at the BMC? Why is the sum Ri + R 2 of the 
gage outputs nearly proportional to M at the BMC? These assertions will be explained more rigorously in 
a later section of the paper. The advantage can also be understood if the two proposed linear combinations 
of the gage outputs are plotted versus the loads at the BMC. Figure 5a shows the difference of the gage 
outputs plotted versus the applied normal force. Figure 5b shows the sum of the gage outputs plotted 
versus the applied moment. In both cases we see that the selected responses are nearly proportional to the 
corresponding load that is applied. Therefore, any regression model developed for the difference and the 
sum of the gage outputs will be dominated by a single linear term if the force and moment at the BMC are 
used as independent variables. This characteristic greatly simplifies the development of a regression model 
of the calibration data. - In the next section of the paper the regression analysis of the difference and the 
sum of the balance gage outputs will be discussed in more detail. 

IV. Regression Analysis of Calibration Data 


A. General Remarks 

The selection of suitable regressors and responses of the sting balance calibration data set has a sig- 
nificant influence on the overall quality of the regression analysis of the data. It was decided to use the 
applied force F and moment M at the BMC as independent variables for the regression analysis. In that 
case, as explained above, (i) the difference R 2 — Ri and (ii) the sum Ri + R 2 of the moment gage out- 
puts should be used as responses. Two different regression model types were selected for each response 
in order to demonstrate the benefit of using an optimized regression model for the analysis of the balance 
calibration data. The first regression model type is a traditional quadratic . The second regression model 
type is the optimized regression model that the candidate search algorithm selected. At first, the regression 
analysis results for the traditional quadratic are discussed. Afterwards, these results will be compared with 
corresponding results that were obtained using the optimized regression models of the calibration data. 

B. Traditional Quadratic 

A traditional quadratic may be selected for the regression analysis of the difference of the gage outputs. 
Then, the difference may be modeled in the least squares sense using the following math model: 


Quadratic 1 : R2 — Ri = 00 + • F + (*2 ■ M + 03 ■ F 2 + 04 ■ M 1 + o 5 ■ F ■ M (la) 


Similarly, the sum of the gage outputs may be modeled using the following quadratic: 


Quadratic 2 : R\ -|- R2 — /?o “I - ft. 1 * F T $2 * AT -|- j3 3 ■ F ^ -|- j3 4 ■ M^ 1 -|- j3§ ■ F ■ M ( 15 ) 


Figure 6a shows regression analysis results for the difference of the gage outputs if the traditional 
quadratic is used for the analysis. This math model has variance inflation factors (VIFs) that are significantly 
larger that the literature recommended liberal threshold of 10 (see list of VIFs printed in red in Fig. 6a). The 
math model described in Eq. (la) has massive near-linear dependencies. Therefore, its predictive capability 
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is put into question even though the standard deviation of the fitted residuals of the responses is only 0.0836% 
of the largest response magnitude. 

Figure 6b shows corresponding regression analysis results for the sum of the gage outputs if a traditional 
quadratic is used for the analysis. Again, the VIFs are large indicating the presence of massive near-linear 
dependencies. Consequently, the predictive capability of the regression model given in Eq. (lb) is also put 
into question even though the standard deviation of the fitted residuals of the responses is only 0.0393% of 
the largest response magnitude. 

C. Optimized Regression Model 

In the next phase of the data analysis the regression model optimization algorithm defined in Fig. 1 
was applied to the difference and the sum of the gage outputs. First, a lower and an upper bound for the 
candidate math model search had to be defined. It was decided to use a simple linear math model as the 
lower bound and a traditional quadratic as the upper bound for the candidate math model search. Then, 
we get for the difference of the gage outputs the following bounds: 

{ Lower Bound 1 : po + Pi ■ F (2a) 

Upper Bound 1 : ao + a\ • F + a<2 • M + <23 ■ F 2 + a<± ■ M 2 + • F • M (2 b) 

The search was performed using the most conservative thresholds for both the primary and secondary 
search constraint: (1) p-value of the ^-statistic threshold => 0.0001 and (2) VIF threshold => 5. A total 
number of 20 regression models were tested during the search. The optimized regression model, i.e., the 
recommended math model, was obtained after the completion of the search. It has the following form: 


Optimized Model 1 : R 2 — Ri = eo + ei • F (2c) 


Figure 7a shows the regression analysis results for the difference of the gage outputs if the optimized 
regression model defined in Eq. (2c) is used for the analysis. The regression model has no near-linear 
dependencies as the VIFs are significantly smaller that the liberal threshold of 10. In addition, only two 
terms are used (instead of six for the traditional quadratic) . The standard deviation of the fitted residuals of 
the difference is only 0.1056% of the largest response magnitude. This value is only slightly larger than the 
value that was obtained for the traditional quadratic (see Fig. 6a). Therefore, it can be concluded that the 
traditional quadratic, i.e., Eq. (la), significantly overfitted the difference of the gage outputs. Overfitting of 
experimental data should be avoided as it can greatly reduce the predictive capability of a regression model 
at data points that were not used for the regression analysis. 

The regression model optimization was applied to the sum of the gage outputs in a similar fasion. Again, 
the lower and upper bound for the candidate math model search were selected to be: 

{ Lower Bound 2 : Ao + Ai ■ M (3a) 

Upper Bound 2 : /3q + Pi • F + @2 * M + Ps ■ F 2 + p± ■ M 2 + /?s ■ F ■ M (3b) 

In the next phase of the analysis the candidate math model search algorithm was applied to the sum of 
the gage outputs. In that case, the algorithm chose the following recommended math model: 


Optimized Model 2 : Ri + R 2 = Po + Pi * M + 772 ■ M 2 ( 3 c) 


Figure 7b shows regression analysis results for the sum of the gage outputs if the optimized regression 
model defined in Eq. (3c) is used for the analysis. The regression model has no near-linear dependencies as 
the VIFs are significantly smaller that the liberal threshold of 10. In addition, only three terms are used 
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(instead of six for the traditional quadratic). The standard deviation of the fitted residuals of the difference 
is only 0.0472% of the largest response magnitude. This value is only slightly larger than the value that 
was obtained for the traditional quadratic (see Fig. 6b). Again, it can be concluded that the traditional 
quadratic, i.e., Eq. (16), significantly overfitted the sum of the gage outputs. 

It is a surprising result of the regression model search that (i) the optimized model of the difference of 
the gage outputs only depends on the intercept and the normal force and that (ii) the optimized model of 
the sum of the gage outputs is only a function of the intercept , the moment , and the square of the moment . 
An explanation of these search results had to be found. 

Fortunately, the overall geometry and design of the sting balance is very simple. Therefore, it will be 
shown analytically in the next section of the paper that the optimized regression model of the difference of the 
gage outputs (i.e., Eq. (2c)) can be obtained after applying classical strain-gage balance load transformations 
to the sting balance. 

It is suspected that the term selection of the optimized regression model of the sum of the gage out- 
puts (i.e., Eq. (3c)) is primarily influenced by the elastic deformation characteristics of the sting balance. 
Therefore, it was decided to use the equations of the deflection of a cantilever beam in order to develop 
an expression for the connection between the balance loads at the BMC and the sum of the gage outputs 
that can be compared with the corresponding optimized regression model. Results of this second theoretical 
analysis will also be presented in one of the next sections of the paper. 


V. Theoretical Regression Model of the Difference of the Gage Outputs 


In this section it is shown how the optimized regression model of the difference of the sting balance 
gage outputs (Eq. (2c)) can be obtained if the classical strain-gage balance load transformation equations 
are applied to the sting balance. At first, an analysis of the connection between the balance loads needs to 
be performed. It is concluded from Fig. 8 that the moment at the BMC is given by the following expression: 


M 


F ■ p 


(4) 


where F is the normal force that is applied in the vicinity of the free end of the sting balance. The sting 
balance has moment gages, i.e., the gage outputs are proportional to moments at the gage location. Therefore, 
in order to develop the relationship between the gage outputs and the force and moment at the BMC, the 
moment at the forward and aft gage needs to be expressed as a function of the loads at the BMC. Using 
again Fig. 8, the moments Mi and M 2 at the center of the forward and aft gage may be written as: 


Mi 


Then, after inserting Eq. (4) into Eqs. (5a) and (56), we get: 
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(5a) 


(56) 


Mi 

F 

= M d 

2 

(6a) 

m 2 

F 

= M+—.d 

(66) 


Equations (6a) and (66) are the classical strain-gage balance load transformations for a moment type 
balance assuming that the BMC is located halfway between the forward and aft gage (see enclosed appendix 
for a completed derivation of the classical strain-gage balance load transformations). Equations (6a) and 
(66) are also a linear system of two equations in two unknowns that may be used to obtain the relationship 
between the loads F and M at the BMC and the moments Mi and M 2 at the electrical centers of the gages. 
After some algebra, we get the solution of this linear system: 


F = 


M2 — M\ 

d 


( 7 ) 
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M 


Mi + M2 
2 


(8) 


As expected, Eqs. (7) and ( 8 ) agree with Eqs. (A. 24a) and (A. 246) for a moment balance that are derived 
in great detail in the appendix of the present paper. Now, the relationship between the loads at the BMC 
and the moments at the gages is known. It remains to find the relationship between (i) the electrical outputs 
and (ii) the moments at the gages. As mentioned above, it is simply assumed that the gage outputs are 
proportional to the moments at the electrical center of the gages. Then, we get the following relationships: 


R 1 

= 71 ■ Mi 

(9a) 

Ri 

= 72 ■ M<2 

(96) 


The sting balance calibration data set may be used to get the numerical values for the two constants of 
proportionality. A least squares analysis of the calibration data leads to the following results: 


7i 

72(^2, M2) 


2.1643 


2.1774 


ftV/V I 
in — lbs 

where 

a = 0.0046 

VV/V I 
in — lbs 

fjiV/V 1 
in — lbs 

where 

cr = 0.0044 

fiv/v 1 

in — lbs 


( 10 a) 

(106) 


Comparing the right hand sides of Eqs. (10a) and (106) we see that the first derivatives of the two gage 
outputs are almost identical. Therefore, the following assumption can be made: 


7 = 7i = 72 

The constant 7 may be computed by fitting a line to the combined calibration data set. We get: 


7 (iJi,i? 2 ,Mi,M 2 ) 


2.1709 


yv/v 

in — lbs 


where a = 0.0079 


yy/y 

in — lbs 


(lla) 


(116) 


Consequently, the final relationship between the electrical output and the moments at the two gages becomes: 


R 1 

= 7 • Mi 

( 12 a) 

R 2 

= 7 ■ m 2 

(126) 


Now, the difference between the two gage outputs can be computed by subtracting the right hand side 
of Eq. (12a) from the right hand side of Eq. (126). Then, we get: 


R 2 — R\ 


7 . 


M^ — Mi 


(13 a) 


The difference of the moments on the right hand side of Eq. (13a) still needs to be replaced. Rearranging 
terms in Eq. (7) we know: 

M 2 - Mi = F-d (136) 

Finally, after replacing the difference of the moments in Eq. (13a) with the right hand side of Eq. (136), 
we get the relationship between (i) the difference of the gage outputs and (ii) the loads at the BMC: 


R 2 — Ri = c 1 m F (14a) 

e'i = 7 • d (146) 
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Now, the right hand side of Eq. (14a) may be compared with the right hand side of Eq. (2c). In both 
equations the applied calibration force F is the first regressor of the sting balance calibration data set. It is 
important to remember that Eq. (14a) was derived assuming that the electrical outputs of the two gages do 
not have a systematic error. Therefore, the intercept term in the theoretical math model, i.e., Eq. (14a), is 
zero. The optimized regression model, i.e., Eq. (2c), was obtained using the actual data set. Therefore, the 
intercept term in the optimized regression model is non-zero as it picked up any systematic error that was 
present in the measured gage outputs. 

VI. Theoretical Regression Model of the Sum of the Gage Outputs 

A. Exact Regression Model 

The equations of the deflection of a cantilever beam under load may be used to show that the optimized 
regression model of the sum of the gage outputs, i.e., Eqs. (3c), can also be derived from first principles of 
engineering mechanics. At first, the sum of the gage outputs needs to be computed. Adding the right hand 
sides of Eq. (12a) and Eq. (126), we get: 

R\ + R% = 7 ‘ Afi + M 2 (15a) 

The sum of the moments on the right hand side of Eq. (15a) still needs to be replaced. Rearranging 
terms in Eq. ( 8 ) we know: 

Mi + M 2 = 2 ■ M (156) 

Now, after replacing the sum of the moments in Eq. (15a) by using the right hand side of Eq. (156), we 
get the relationship between the sum of the gage outputs and the loads at the BMC: 

Ri + R 2 = [ 7 • 2 1 • M (15 c) 

The moment at the BMC equals the normal force times the moment arm. Therefore, the moment at 
the BMC can be expressed as follows: 

M = F ■ Z (16) 

The moment arm Z is the true moment arm of the normal force relative to the BMC assuming that the 
sting balance experiences an elastic deformation. This moment arm can be expressed as the sum of (i) a 
“nominal” moment arm Z 0 at zero force (equals the moment arm for a rigid sting) and (ii) a moment arm 
correction AZ that takes the deformation of the sting balance under load into account. Then, we get: 

Z = Z 0 + AZ (17) 

Inserting the right hand side of Eq. (17) into Eq. (16), we get: 

M = F ■ [ Zo + AZ 1 (18) 

Now, the moment M in Eq. (15c) is replaced using the right hand side of Eq. (18). Then, Eq. (15c) becomes: 

Rx + R 2 = [ 7 ■ 2 1 • F • [ Zo + AZ 1 = [ 7 • 2 1 ■ F ■ Z G + [ 7 ■ 2 1 ■ F ■ AZ (19) 

During the balance calibration the “nominal” moment at the BMC was determined by computing the 
product of the applied calibration force and the “nominal” moment arm. The “nominal” moment is the 
second regressor of the calibration data set. It can be expressed as follows: 

M 0 = F-Z 0 (20) 

Now, after using Eq. (20) in order to introduce the “nominal” moment in Eq. (19), we get: 
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i?l + i?2 


7-2 


■M 0 + 


7-2 


• F • Al 


( 21 ) 


The change of the moment arm A l due to the elastic deflection of the sting balance under load still 
needs to be determined. It is assumed that the moment arm change is primarily caused by the bending 
of the sting balance under load. The sting balance itself may be considered as a cantilever beam in order 
to quantitatively assess the moment arm change due to bending. The exact solution of this moment arm 
correction due to sting deformation is derived in the next section of the paper. 

B. Exact Solution of Moment Arm Correction 

Figure 8 shows the loads that act on the sting balance during calibration. A calibration fixture is 
attached to the free end of the sting. The fixture has two tasks: (i) it allows for the application of the 
normal force at a desired “nominal” moment arm and (ii) it makes it possible to attach an AMS unit to the 
balance that is needed for the leveling of the calibration fixture after the load is applied. It is important to 
remember that the calibration force is applied at locations on the calibration fixture that do not necessarily 
coincide with the free end of the sting (green dot in Fig. 8 marks the location of the calibration force on the 
calibration fixture). Therefore, in the general case, two types of beam bending need to be superimposed in 
order to estimate the moment arm change. The first bending is caused by a concentrated force (see Fig. 9a) . 
The second bending is caused by a moment at the free end of the sting (see Fig. 9b) . Consequently, the total 
moment arm correction is the sum of the superimposed corrections due to the two types of bending. We get: 

A l = l\ + (2 2) 

It remains to find a relationship between the beam deflections 5\ and 62 and the moment arm corrections 
Zi and Z 2 . Similar triangles in Fig. 9a and Fig. 9b may be used for this purpose. We get: 

l\ _ 1 2 _ h 

Si 62 l 0 /cos(j) 

Now, after solving Eq. (23) for the moment arm correction Zi , we get: 

^ l 0 /cos(j) ^ 

It is suspected that the deflection 5\ due to the concentrated force is significantly larger than the 
deflection 62 due to the moment at the free end. The equations of maximum cantilever beam deflection 
given on p.170 of Ref. [8] may be used to investigate this hypothesis. We know, using these equations, that 
the ratio between 5\ and 82 may be approximated as follows: 


(24) 


(23) 


Si 

[ F ■ coscj) ] ■ [l 0 /cos(p] 3 


2 El 

2 

lo 

62 ~ 

3 El 


[s • F ■ sin(j) ] ■ \l 0 /cos(f)\ 

~3~ 

s ■ sin (f) 


^ ^ ^ Vi ^ s 

concentrated force moment at free end 


The estimated bend angle <ft of the sting in Fig. 9a and Fig. 9b is between 15° and 20°. The moment 
arm Z D is on the order of 20[in]. In addition, the distance s in Fig. 9b is small (^ 2 [in]). Therefore, we get 
for the ratio of the maximum deflections the estimate: 

20 < A < 26 (256) 

02 

Consequently, we conclude that the deflection Si is significantly larger than the deflection £ 2 . This 
conclusion also applies to the moment arm corrections. We get: 

|^i| > | <fe| =>• |ti | > \h\ => h ~ 0 (26) 
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Finally, after (i) replacing l\ and h in Eq. (22) using Eq. (24) and the result of Eq. (26), and after 
(ii) using the symbol 6 instead of , we get for the total moment arm correction the equation: 


A l 


h 

lo/coscj) 


■5 


(27) 


An estimate of the beam deflection 6 due to a concentrated force is still needed so that the coefficients 
of the theoretical regression model of the sum of the gage outputs can be quantified. The deflection can be 
found using the Area-Moment Method (see Ref. [8] for more detail). Then, the deflection equals 


<5 


r b m{x) ■ x 

« E(X) ■ I(X) 


(28) 


where x is the beam coordinate, a and b are the coordinate of the fixed and free end of the beam, M (x) is 
the local moment, E(x) is the modulus of elasticity, and I(x) is the moment of inertia of the cross-section 
of the beam. The modulus of elasticity and the moment of inertia need to be moved in front of the integral 
symbol so that the deflection can be estimated. This can be done by using a generalized form of the Mean 
Value Theorem of the Integral Calculus (see Ref. [9], pp. 126-128). The theorem can be described as follows: 


If f(x) and p(x) are continuous functions in a < x < b, and p(x) > 0 or, 
more general, p(x) does not change sign in the interval a < x <b, then 

[ f(x)-p(x)dx = /(£)• t p(x) dx (29) 

J a J a 

where £ is some constant that fulfills the condition a < £ <b. 


The following four substitutions have to be made so that Eq. (29) may be applied to Eq. (28): 


a = 0 (30a) 

b = l Q /cos(f) (306) 

= ew- m {m 

p(x) = M(x)-x (30 d) 

Then, Eq. (28) can be written in the following form: 


6 = 


L 


l Q /cOS(f> 


m) ■ m 

From Fig. 9a we conclude that the following equation is valid: 


M(x) • x dx 


M(x) ■ x = 


F • cos 4> • x 


■ x 


(31) 


(32) 


In the next step, after replacing the integrand in Eq. (31) by the right hand side of Eq. (32) and knowing 
that F ■ cos 0 is a constant concentrated load at the assumed free end of a cantilever beam, we get: 


6 = 


F ■ cos 4> 

m ■ m 


L 


lo/cOS(j> 


x 2 dx 


(33a) 
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Finally, after solving the integral in Eq. (33a) analytically, we get the following value for the deflection 
of the sting balance: 

6 = 3 ■ E{£) ■ J(0 ■ cos 2 0 

Now, it is possible to get the moment arm correction due to bending. It is only required (i) to replace 
the deflection 6 in Eq. (27) with the right hand side of Eq. (336). Then, the moment arm correction becomes: 


A l = 


h 


6 = 


h 


F-l 3 


lo/coscp l 0 /cos<j> 3 ■ E(£) ■ I(£) ■ cos 2 (f) 

Equation (34) can be simplified further. After some algebra we get for the moment arm correction: 


(34) 



C. Approximated Solution of Moment Arm Correction 

Unfortunately, the numerical value of the moment arm correction due to bending cannot easily be 
computed using Eq. (35) as the exact value of E(£) • /(£) is unknown. However, it is possible to develop 
an estimate for E(£) • /(£) by using an alternate expression for the deflection of a cantilever beam. From 
Ref. [8], p.159, we also know that 

dQ _ M(x) 


dx E( x ) • I(x) 

where 0 is the deflection angle. Then, Eq. (28) can also be written as follows: 


/ 


M(x) • x 
E(x) ■ I(x) 


dx 


f 


dQ 

dx 


• x dx 


(36) 


(37) 


In most cases the slope dQ/dx is a function of the beam coordinate x. A simplifying assumption, 
however, can be made that will make it possible to assess the magnitude of the coefficients of the regression 
model of the sum of the gage output. It is simply assumed that the slope dQ/dx is constant: 


dQ 


constant 


0(5) -0(a) 


dx b — a 

Then, the deflection of the cantilever beam at the free end can be approximated as follows: 


0(5) - 0(q) 
5 — a 


■f 


x dx = 


0(5) -0(a) 5 2 — a 2 


0(5) -0(a) 


a T 5 


5 — a 2 

Equation (39) may be applied to the sting balance if the following assumptions are made: 
a = 0 ; 5 = l 0 /cos(j) ; 0(a) = 0 ; 0(5) ==> from experiment 
Then, the deflection at the free end becomes: 

Now, after replacing 6 in Eq. (336) with the right hand side of Eq. (40), Eq. (336) becomes: 

Q{l 0 /cos4>) F ■ l 3 


(38) 


(39) 


lo/coscp ■ 


3 ■ E(£) ■ /(£) ■ cos 2 4> 


(40) 


( 41 ) 


12 

American Institute of Aeronautics and Astronautics 




After some algebra we get for the inverse of the product of the modulus of elasticity with the moment 
of inertia the following approximation: 


1 

W) • m 


3 • Q(l 0 /cos(f)) 

2 • F • Z 0 • l 0 /cos(j) 


From Fig. 8 we also know that: 

lo/coscf) = \J h 2 + 1*1 


(42) 

(43) 


Then, using Eq. (20) to replace F ■ l 0 in Eq. (42) and after replacing l 0 jcosf> with the right hand side 
of Eq. (43), Eq. (42) can be written as: 


m-m 


3 • Q(l 0 /cos(f)) 


(44) 


Fortunately, the deflection angle at the free end as a function of the applied loads was determined during 
the preparation of the sting balance for a test at Ames Research Center. This data was fitted using a least 
squares fit and the following regression model of the deflection angle at the free end was obtained: 


Q(l 0 /cos(j)) 


ki - M 0 + k 2 * F 


180° 


from experiment 


where 


h 


k 2 


+0.00329615 


-0.01702055 


deg 


in — lbs 
deg 


lbs 


The force F in Eq. (45a) can be replaced by using the “nominal” moment and moment arm. We get: 

F = Mo 

Then, Eq. (45a) can be written in the following form: 


lo 


Q(lo/cOS(f)) 


fci • M 0 + (k 2 /lo) * Af 0 


180° 


+ ( k 2 /lo ) 


180° 


• M 0 


(45 a) 

(456) 

(45c) 

(46) 

(47) 


Equation (47) may be simplified further. The range of the “nominal” moment arm is known from 
balance calibration records. Therefore, the following relationship applies: 


16.52 [in] < lo < 25.92 [in] 

Now, using the result of Eq. (48) with Eq. (45c), we get the following estimate: 


0.00066 


deg 


in — lbs 


< \k 2 /lo\ < 0.00103 


deg 


in — lbs 


(48) 


(49) 


Comparing the absolute value of the coefficient k\ given in Eq. (456) with the range of the absolute 
value of the term k 2 /l 0 we see that \ki \ is three to five times larger than the term \k 2 /l 0 \. Therefore, we can 
make the following simplifying assumption: 


\ki\ » \k 2 /l 0 \ => 

Then, using the result of Eq. (50), Eq. (47) becomes: 

Q(l 0 /cos4>) ki 


\k 2 /l 0 


0 


180° 


■M 0 


(50) 


( 51 ) 
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Finally, after using the result of Eq. (51) in Eq. (44), we get: 


ki ■ 7 r 


E (0 ■ HO 120° ■ V h? + 

Now, using the result of Eq. (52) in Eq. (35), we get for the moment arm correction the equation: 

ki-TT-h ( 1 \ F /2 


A l » - , 

360° • sj h? + V cos <P 

We also know from Fig. 8 that the following relationship applies: 

h 


sin 4> = 


Therefore, Eq. (53) can be expressed as follows: 


V v + % 


(52) 


(53) 


(54) 



D. Approximated Solution of Exact Regression Model 

Finally, the moment arm correction estimate defined in Eq. (55) is inserted into Eq. (21). Then, the 
sum of the gage outputs becomes: 


Ri + R<2 


7 * 2 


■M 0 + 


7' fci ' 
180° 


■ tan (j) 


F 2 -It 


(56) 


The numerical value of tancp in Eq. (56) depends on the “nominal” moment arm l Q . We know that 

tan 0 => tan 0(/ o ) (57) 

We can use again the definition of the “nominal” moment (Eq. (20)). Then, after substituting F ■ Z D in 
Eq. (56) using M 0 , we get the final form of a theoretical regression model of the sum of the gage outputs: 



In the next part of the paper the theoretical regression coefficients of the regression models for the 
difference and the sum of the gage outputs are compared with results of the optimized regression models. 

VII. Regression Coefficient Comparison 


A. General Remarks 

In the previous sections exact and approximate solutions of the regression models of (i) the difference 
and (ii) the sum of the gage outputs of the sting balance were developed. The theoretical and fitted values 
of the regression model coefficients need to be compared in order (i) to assess the validity of the assumptions 
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that were used to develop the theoretical regression models and (ii) to arrive at a better understanding of the 
calibration data. The theoretical values of the regression coefficients are computed using balance geometry 
data and other information. The regression cofficients of the optimized regression models, on the other hand, 
are taken directly from the listed regression analysis results (column three in Figs. 7a and 7b). 

B. Difference of Gage Outputs 

At first, the optimized regression model of the difference of the gage outputs is investigated. The 
optimized regression model is defined in Eq. (2c). The parameter d, i.e., the distance between the two 
moment gages of the sting balance is needed in order to determine the theoretical value of the first coefficient 
ei that is related to the applied calibration force F. The theoretical value of the coefficient is defined in 
Eq. (146). From engineering drawings of the sting balance we know that the parameter d is given as: 

d = 5.5 [in] (59) 

Then, using the estimate for 7 given in Eq. (116) and the value of the distance d given in Eq. (59), the 
theoretical value of the coefficient becomes: 



The corresponding coefficient value for the optimized regression model is listed in the table in Fig. 7a 
(column three). The coefficient has the following value: 



Comparing the right hand side of Eq. (60a) with the right hand side of Eq. (606) we see that the 
theoretical and the fitted value show excellent agreement. 

C. Sum of Gage Outputs 

In the next step the coefficients of the optimized regression model of the sum of the gage outputs needs 
to be compared. The optimized regression model is defined in Eq. (3c). At first, the coefficient of the linear 
term is investigated. Its theoretical value is given in Eq. (586). It has the following solution: 



The corresponding coefficient value for the optimized regression model can be found in the third column 
in Fig. 7b. It is given as: 



Again, after comparing the right hand side of Eq. (61a) with the right hand side of Eq. (616), we see 
that the theoretical and the fitted values show excellent agreement. 

It remains to compare the coefficient of the square term of the regression model of the sum of the 
gage output. The estimated theoretical value of the term is given on the right hand side of Eq. (58c). 
Unfortunately, the angle 0 is not constant for the calibration data set. Therefore, the range of the magnitude 
of the theoretical value of the coefficient needs to be computed as a function of the range of angle 0. From 
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the calibration records we know that the range of angle 0 is given as follows: 


15° < 0(Z O ) < 20 ° (62a) 

Therefore, tan 0 has the following range: 

0.27 < tan <j>(l 0 ) < 0.36 (626) 

Consequently, after using (i) the range given in Eq. (626) and (ii) the values for 7 and k\ given in 
Eq. (116) and Eq. (456) in Eq. (58c), we get the following range for the theoretical value of the second 
coefficient of the regression model of the sum of the gage outputs: 


Theoretical Analysis : 3.4 x 10 5 

(#*V/V) j 

(in — lbs) 2 _ 

< 772 < 4.5 x 10 “ 5 

W/v) 1 

_(in — lbs) 2 

(63a) 


Again, the corresponding value of the second coefficient of the optimized regression model can be found 
in the third column in Fig. 7b. It is given as: 


Least Squares Fit : 772 


-7.3 x 1(T 6 


w 0 

(in — lbs) 2 


(636) 


Two observations can be made if the range of the theoretical value (Eq. (63a)) is compared with the 
least squares estimate (Eq. (636)): (1) the absolute value of the theoretical value is approximately four to 
six times larger than the fitted value; (2) the sign of the theoretical value does not equal the sign of the 
fitted value. How can these differences between the theoretical estimate and the fitted value of the second 
coefficient be explained? It must be remembered that the calculation of the theoretical estimate did not take 
the fact into account that the calibration fixture was leveled each time after a specific calibration load was 
applied. Figure 10 shows the impact of the leveling of the calibration fixture on the moment arm. Several 
observations can be made using Fig. 10: 

Observation 1: The leveling of the balance partially counteracts the increase of the moment arm due to 
the elastic deformation of the sting balance under load. Therefore, the magnitude of the fitted coefficient 
must be significantly smaller than the magnitude of a theoretical estimate of the coefficient as the theoretical 
estimate does not account for the leveling of the calibration fixture. 

Observation 2: The final position of the force attachment point (i.e., position 3 in Fig. 10) does no longer 
coincide with the original position of the load attachment point (i.e., position 1 in Fig. 10) because the sting 
balance remains elastically deformed even though the calibration fixture was leveled. 

Observation 3: The unloaded sting is defined using the line ABC in Fig. 10. This line was also used to 
define the “nominal” moment arm Z Q that is used to compute the “nominal” moment M 0 . The loaded sting 
(after leveling) is defined using the line A f B f C. The leveling removed most of the moment arm change due 
to sting bending. In addition, due to the remaining elastic deformation of the sting under load, the sting 
segments A f B f and B f C are no longer straight lines. Therefore, the total length of the moment arm must 
have actually been reduced after the calibration fixture was leveled under load. This explains the observation 
that the sign of the fitted second coefficient of the regression model is negative and not positive. The sign 
of the coefficient would have been positive if the leveling of the calibration fixture would not have removed 
the moment arm change due to sting bending. 

VIII. Summary and Conclusions 

A regression model optimization algorithm was applied to wind tunnel sting balance calibration data. 
The optimization algorithm correctly predicted that (i) the difference of the gage outputs should be modeled 
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using an intercept term and the normal force at the BMC and that (ii) the sum of the gage outputs should 
be modeled using an intercept term, the moment at the BMC, and the square of the moment at the BMC. It 
was also shown that the difference and the sum of the gage outputs best describe the physical characteristics 
of a sting balance if the normal force and the moment at the BMC are the chosen calibration loads. 

Classical strain-gage balance load transformations were applied in order to demonstrate that the opti- 
mized regression model of the difference of the gage outputs is supported by a rigorous analysis of the physics 
of the balance. The theoretical analysis of the relationship between the loads at the BMC and the sum of 
the gage outputs was more difficult. In that case, the equations of the deflection of a cantilever beam had to 
be applied in order to confirm the validity of the optimized regression model. It was also demonstrated that 
the leveling of the sting balance during the calibration is important. It removes most of the moment arm 
change that is introduced in the calibration data set due to the bending of the sting balance under load. The 
analysis also showed that the remaining elastic deformation of the leveled sting balance is mostly responsible 
for the presence of the square of the moment in the regression model of the sum of the gage outputs. 

In conclusion, the regression analysis of the sting balance calibration data set is a rare example of a 
situation when the regression model of a strain-gage balance calibration data set can be derived from first 
principles of physics and engineering. It is also interesting to note that the regression model optimization 
algorithm predicted the same math term combination for the balance calibration data using only a set of 
statistical quality metrics. 
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Appendix: Classical Strain-Gage Balance Load Transformations 


A. Introduction 

Different types of strain-gage balance designs are used in wind tunnel testing. A force type balance, 
for example, is designed such that a balance gage responds to a force at the gage location. A moment 
type balance, on the other hand, has gages that respond to a moment at the gage location. An alternate 
derivation of the classical load transformation equations for these two types of balances is presented in this 
appendix. These transformation equations describe the relationship between loads at the balance moment 
center (BMC) and loads at the location of the gages. 

B. Force Type Balance 

Figure 11 shows forces and moments that act on a typical strain-gage balance. The load transformation 
equations for a force type balance can easily be derived if the force F and moment M at the BMC are 
replaced by two forces F\ and F 2 at the location of the gages. The sum of these two forces equals the force 
at the BMC. We get: 


F = F 1 + F 2 (A.l) 

Similarly, we know that the total moment at the BMC equals the sum of the moment contributions 
from the two forces that replace the force and moment at the BMC. It is assumed that the location of the 
BMC and the location of the gages are described in a coordinate system that uses the balance centerline as 
the coordinate system axis (see Fig. 11). Then, the BMC has the coordinate /i, the center of the forward 
gage the coordinate Ci, and the center of the aft gage the coordinate C 2 . Now, two coordinate differences d\ 
and d 2 may be defined that describe the distance between the BMC and the forces F\ and F 2 at the balance 
gage locations. We get: 


d\ = fL — ci (A. 2a) 

d 2 = fL - c 2 (A. 2b) 

Then, using the sign definitions of (i) the balance loads and of (ii) the coordinates that are depicted in 
Fig. 11, the total moment at the BMC can be expressed as follows: 

M = Fi - di - F 2 - (-d 2 ) (A. 3) 


From Eq. (Al) we know that: 


F 2 = F - F 1 (AA) 

The right hand side of Eq. (A A) may be used to replace the force F 2 in Eq. (A. 3). We get: 


M = F 1 • di - (F - Fi) - (-da) 
Rearranging terms in Eq. (A. 5) and after some algebra we get: 


Fi 


(-da) 

di — d 2 


+ 



Similarly, using Eq. (A. 6) to replace force F\ in Eq. (AA) and after some algebra, we get: 


(A5) 


(A6) 


F 2 


F ■ 


d i 

d\ — d 2 



(A7) 


Finally, using Eqs. (A. 2a) and (A. 2b) to replace the moment arms in Eqs. (A. 3), (A. 6), and (A. 7), we get 
the set of load transformation equations for a force type balance: 
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C. Moment Type Balance 

A moment type balance is designed such that a balance gage responds to a moment at the gage location. 
Therefore, we need to express the moment at the gage location as a function of the force and moment at the 
BMC in order to derive the desired load transformation equations. We know from Fig. 12 that the moment 
Mi at the forward gage is caused by the force at the aft gage. Therefore, using the gage coordinates 
defined in Fig. 12, we can write: 

Mi = (—1) • F 2 • (c 2 — ci) 

Now, using the right hand side of Eq. (A.8d) to replace the force F 2 in Eq. (A. 9), we get 

M 1 = (-1) ■ ( F • M ~ Cl - M 1 V (C 2 - ci) 

V 02 - Cl c 2 - Cl ) 

Finally, after simplifying and rearranging the right hand side of Eq. (A. 10), we get: 

Mi = M - F-(fi - c i) (A. 11) 

The derivation of the moment M 2 at the aft gage uses the same approach. This time, the moment M 2 at 
the aft gage is caused by the force F\ at the forward gage. We get: 

M 2 = (+1) • Fi • (c 2 — ci) 

Now, using the right hand side of Eq. (A. 8c) to replace the force Fi in Eq. (A. 12), we get 

M 2 = (+1) • ( F • C2 ~ M + M V (c 2 - Cl) 

V 02 - Cl c 2 - Cl ) 

Finally, after simplifying and rearranging the right hand side of Eq. (A. 13), we get: 

M 2 = M + F-(c 2 — //) (A. 14) 

The remaining transformations can be obtained from Eqs. (A. 11) and (A. 14). Subtracting the left and right 
hand side of Eq. (A. 11) from the left and right hand side of Eq. (A. 14), we get: 

M 2 — Mi = F • (c 2 — ci) 

Then, after solving Eq. (A. 15) for the force at the BMC, we get the following relationship 

M 2 — Mi 
F - 0, - ft 
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(A. 15) 
(A. 16) 


(A. 12) 
(A.13) 


(A9) 

(A10) 




In addition, after solving Eq. (All) for the moment at the BMC, we also know that: 


M = Mi + F • ({i — ci) 

Then, using the right hand side of Eq. (A 16) to replace the force F in Eq. (A 17), we get: 


M 


Mi + 


/ M 2 - Mi \ 

V 02 - Cl ) 


■ (M - ci) 


We also know that the right hand side of Eq. (A 18) can be written in the following form: 


(A 17) 


(A 18) 


Ml +(— ? — Y(/x - ci) = Mr f— — )+M 2 -(— —) - Mi- f— — ) (A19) 

\ c 2 - Cl J V C 2 - Cl ) \ C 2 - Cl ) \C 2 - Cl ) 

Now, after (i) simplifying the right hand side of Eq. (A 19) and after (ii) replacing the right hand side of 

Eq. (A18) with the result, we get for the moment M at the BMC: 

M = Afi • ( - - M ) + M 2 • ( - - - ) (A20) 

\ C 2 ~ Ci ) V °2 - Cl J 

Finally, the set of load transformation equations for a moment type balance can be summarized: 


II. LOAD TRANSFORMATIONS - MOMENT TYPE BALANCE 

(univeral relationships - valid for any given balance moment center location) 

TP — 

to 

1 

£ 

(A21a) 

jp — 

C*2 ~ Cl 

M = 

Ml .( C2 -M +m 2 .( 

V c 2 - ci y \c 2 - ci j 

(A216) 

Mi = 

M — F • (fi — ci) 

(A21c) 

M 2 = 

M + F ■ (c 2 - /i) 

(A21d) 


D. Load Transformation Equations for Simplified Gage Configuration 

The BMC is often located halfway between the forward and aft gage of a balance. In that case, the load 
transformation equations can be simplified significantly if the gage distance d is introduced as a parameter. 
Then, after placing the coordinate system origin at the location of the forward gage, we get: 


Cl 

= 0 

(/ 1.22a) 

C2 

= d 

{A. 22b) 


= d/2 

{A. 22c) 


Now, after applying Eqs. (A 22a), (A226), (A22c) to Eqs. (A8a), (A8&), (A8c), (A8d) and after some 
algebra, we get the simplified transformations for the force type balance: 
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III. SIMPLIFIED LOAD TRANSFORMATIONS - FORCE TYPE BALANCE 

(valid if the balance moment center is located halfway between the forward and aft gage) 


F 

= Fi + F% 

(A. 23a) 

M 

^3 

£ 

1 

GT 

11 

{A. 236) 

F i 

F M 

~2 + ~d~ 

(A. 23c) 

Fi 

F M 

~ ~2 T 

(A. 23d) 


Similarly, after applying Eqs. (A22a), (A226), (A 22c) to Eqs. (A 21a), (A 216), (A 21c), (A21d) and after 
some algebra, we get the simplified transformations for the moment type balance: 


IV. SIMPLIFIED LOAD TRANSFORMATIONS - MOMENT TYPE BALANCE 

(valid if the balance moment center is located halfway between the forward and aft gage) 


TP — 

to 

1 

Ss 

(A. 24a) 

r — 

d 

M = 

Mi + M2 

(A. 246) 

2 

Mi = 

F 

M d 

2 

(A. 24c) 

M 2 = 

F 

M + --d 

(A.24d) 
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Fig. 1 Key elements of candidate math model search algorithm. 
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Fig. 2 Forces and moments acting at different locations on a wind tunnel sting balance. 
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Fig. 3a Sideview of sting balance with attached AMS (angle measurement system) unit. 
(Courtesy of The Gulfstream Aerospace Corporation, Savannah, Georgia.) 



Fig. 3b Location of forward and aft gage of sting balance after removal of gage fairings. 
(Courtesy of The Gulfstream Aerospace Corporation, Savannah, Georgia.) 



Fig. 3c Calibration of sting balance in upside-down orientation using weights. 
(Courtesy of The Gulfstream Aerospace Corporation, Savannah, Georgia.) 
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Fig. 4 Calibration forces and moments at the balance moment center for the four selected moment arms. 



Fig. 5a Difference of gage outputs plotted versus the applied force at the balance moment center. 
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Fig. 5b Sum of gage outputs plotted versus the applied moment at the balance moment center. 
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REGRESSION COEFFICIENT ESTIMATES AND STATISTICAL METRICS (R2-R1) 

REGRESSION MODEL HIERARCHY CHARACTERISTICS - HIERARCHICAL 


MATH TERM 

MATH TERM 

COEFFICIENT 

STANDARD T 

-STATISTIC OF 

P -VALUE OF 

VI F 

VI F 

INDEX 

NAME 

VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 

INTERCEPT 

-0.1051 

+0.0259 

-4.0605 

- 

- 

- 

2 

F 

+ 12.0427 

+0.0196 

+ 613.3586 

< 0.0001 

+ 107.0591 

+ 107.0591 

3 

M 

+0.0083 

+0.0010 

+8.4039 

< 0.0001 

+ 107.0591 

+ 107.0591 

6 

F*F 

+0.0050 

+0.0134 

+0.3704 

+0.7117 

+6558.3359 

+6358.3359 

7 

M*M 

+5.6062e— 06 

+3.1 574e— 05 

+0.1776 

+0.8594 

+5290.8719 

+5290.8719 

10 

F*M 

-0.0004 

+0.0013 

-0.2820 

+0.7784 

+22933.7747 

+22933.7747 

RESIDUAL(R2-R1), % 

OF LARGEST RESPONSE MAGNITUDE 

; STANDARD 

DEVIATION = 0.2044 [microV/V] 

OR 0.0836 


ARITHMETIC MEAN = 

1.02e— 11 % ; RATIO BETWEEN 

LARGEST RESIDUAL AND STANDARD DEVIATION = 

= 4.07 ; LARGEST RESPONSE 

MAGNITUDE = 244.53 [microV/V] 


1.00 % 


0.00 % 


* - 1.00 % 




i g n n n i 

a h n g — — f 



Is 

y — 4J- 


- 100.00 % 


-50.00 % 0.00 % 50.00 % 

R2-R1 (RESPONSE, % OF LARGEST MAGNITUDE) 


100.00 % 


Fig. 6a Gage Output Difference: Least squares solution using a traditional quadratic as the regression model. 


REGRESSION COEFFICIENT ESTIMATES AND STATISTICAL METRICS (R1+R2) 

REGRESSION MODEL HIERARCHY CHARACTERISTICS - HIERARCHICAL 


MATH TERM 

MATH TERM 

COEFFICIENT 

STANDARD 

T— STATISTIC OF 

P -VALUE OF 

VI F 

VI F 

INDEX 

NAME 

VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 

INTERCEPT 

+0.1407 

+0.0889 

+ 1.5827 

- 

- 

- 

2 

F 

+0.31 10 

+0.0675 

+4.6085 

< 0.0001 

+ 107.0591 

+ 107.0591 

3 

M 

+4.3289 

+0.0034 

+ 1270.3971 

< 0.0001 

+ 107.0591 

+ 107.0591 

6 

F»F 

+0.1096 

+0.0459 

+2.3852 

+0.0186 

+6358.3359 

+6358.3359 

7 

M*M 

+0.0002 

+0.0001 

+ 1.5687 

+0.1193 

+5290.8719 

+5290.8719 

10 

F»M 

-0.0091 

+0.0045 

-2.0298 

+0.0445 

+22933.7747 

+22933.7747 


RESIDUAL(R1 + R2), % OF LARGEST RESPONSE MAGNITUDE ; STANDARD DEVIATION = 0.7026 [microV/V] OR 0.0393 % 


ARITHMETIC MEAN = -2.37e-14 % ; RATIO BETWEEN LARGEST RESIDUAL AND STANDARD DEVIATION = 2.57 ; LARGEST RESPONSE MAGNITUDE = 1785.69 [microV/V] 



-100.00 % -50.00 % 0.00 % 50.00 % 100.00 % 

R1+R2 (RESPONSE, % OF LARGEST MAGNITUDE) 


Fig. 6b Gage Output Sum: Least squares solution using a traditional quadratic as the regression model. 
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REGRESSION COEFFICIENT ESTIMATES AND STATISTICAL METRICS (R2-R1) 

REGRESSION MODEL HIERARCHY CHARACTERISTICS = HIERARCHICAL 


MATH TERM 

MATH TERM 

COEFFICIENT 

STANDARD 

T- STATISTIC OF 

P-VALUE OF 

VI F 

VI F 

INDEX 

NAME 

VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 

INTERCEPT 

-0.1116 

+0.0229 

-4.8676 

- 

- 

- 

2 

F 

+ 12.2069 

+0.0024 

+5174.2210 

< 0.0001 

+ 1.0000 

+ 1.0000 


RESIDUAL(R2 — R1), % OF LARGEST RESPONSE MAGNITUDE ; STANDARD DEVIATION = 0.2583 [microV/V] OR 0.1056 % 

^ 1.00 % 

5 

i 

CM 
C£ 

LT 0.00 % 

< 

ZD 

q 
in 
l±j 

K - 1.00 % 

-100.00 % -50.00 % 0.00 % 50.00 % 100.00 % 

R2-R1 (RESPONSE, % OF LARGEST MAGNITUDE) 


ARITHMETIC MEAN = 5.38e-15 % ; RATIO BETWEEN LARGEST RESIDUAL AND STANDARD DEVIATION = 2.91 ; LARGEST RESPONSE MAGNITUDE = 244.53 [microV/V] 


4- h- a r 

^ D n 


1 a B — -H 

i — e — H - ti 


Nis* U i 

L- H--#- H 


Fig. 7a Gage Output Difference: Least squares solution using BALFIT’s optimized regression model. 


REGRESSION COEFFICIENT ESTIMATES AND STATISTICAL METRICS (R1+R2) 

REGRESSION MODEL HIERARCHY CHARACTERISTICS - HIERARCHICAL 


MATH TERM 

MATH TERM 

COEFFICIENT 

STANDARD 

T— STATISTIC OF 

P-VALUE OF 

VI F 

VI F 

INDEX 

NAME 

VALUE 

ERROR 

COEFFICIENT 

COEFFICIENT 

(PRIMARY) 

(ALTERNATE) 

1 

INTERCEPT 

+0.2028 

+0.0996 

+2.0360 

- 

- 

- 

3 

M 

+4.3446 

+0.0004 

+ 1 1 128.4463 

< 0.0001 

+ 1.0000 

+ 1.0000 

7 

M*M 

-7.3372e— 06 

+ 1.7686e— 06 

-4.1486 

< 0.0001 

+ 1.0000 

+ 1.0000 


RESIDUAL(R 1 + R2), % OF LARGEST RESPONSE MAGNITUDE ; STANDARD DEVIATION = 0.8430 [microV/V] OR 0.0472 % 


ARITHMETIC MEAN = 9.39e-14 % ; RATIO BETWEEN LARGEST RESIDUAL AND STANDARD DEVIATION = 2.07 ; LARGEST RESPONSE MAGNITUDE = 1785.69 [microV/V] 



-100.00 % -50.00 % 0.00 % 50.00 % 100.00 % 

R1+R2 (RESPONSE, % OE LARGEST MAGNITUDE) 


Fig. 7b Gage Output Sum: Least squares solution using BALFIT’s optimized regression model. 
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F 



Fig. 8 Forces and moments acting on a sting balance during calibration. 



Fig. 9a Moment arm correction due to force induced beam deflection. 
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CORRECTION 



Fig. 9b Moment arm correction due to moment induced beam deflection. 



Fig. 10 Influence of (i) elastic sting deflection and (ii) calibration fixture leveling on the moment arm. 
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FORWARD 

GAGE 


+F 



Fig. 11 Forces and moments acting at (i) the balance moment center and (ii) the gages of a strain-gage balance. 


+F 



Fig. 12 Moments acting at the forward and aft gage of a strain-gage balance. 
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